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ABSTRACT 

O I Aims. In this paper we confirm the presence of a globally-ordered, kG-strength magnetic field in the photosphere of the young O star 
9 l Orionis C, and examine the properties of its optical line profile variations. 

Methods. A new series of high-resolution MuSiCoS Stokes V and / spectra has been acquired which samples approximately uniformly the 
rotational cycle of 8 1 Orionis C. Using the Least-Squares Deconvolution (LSD) multiline technique, we have succeeded in detecting variable 
^ Stokes V Zeeman signatures associated with the LSD mean line profile. These signatures have been modeled to determine the magnetic field 
• • ■ geometry. We have furthermore examined the profile variations of lines formed in both the wind and photosphere using dynamic spectra. 
m £^ ' Results. Based on spectrum synthesis fitting of the LSD profiles, we determine that the polar strength of the magnetic dipole component 
, is 1150 5> S d S 1800 G and that the magnetic obliquity is 27° S B S 68°, assuming i = 45 ± 20°. The best-fit values for i = 45° are 
$_j ■ Bi = 1300 ± 150 (l<r) G and B = 50° ± 6° (lcr). Our data confirm the previous detection of a magnetic field in this star, and furthermore 
' demonstrate the sinusoidal variability of the longitudinal field and accurately determine the phases and intensities of the magnetic extrema. 
The analysis of "photospheric" and "wind" line profile variations supports previous reports of the optical spectroscopic characteristics, and 
provides evidence for infall of material within the magnetic equatorial plane. 
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1. Introduction 

The detection of magnetic fields in O-type stars has proven 
to be a re markably challenging observational problem (e.g., 
lDonatietal.ll200lh . The apparent absence of magnetic signa- 
tures has often been interpreted as a sel ection effect, s ince 
about 5% of B- and A-type stars (see, e.g., Johnson 2005) do 
display organized magnetic fields with disk-averaged strengths 
rangin g from a few hundred G to several tens of kG ( Mathvs 
I2001I) . These fields are believed to be the fossil remnants of 
either interstellar fields swept up during the star formation pro- 
cess or fields produced by a pre-main sequence envelope dy- 

Send offprint requests to: G.A. Wade, e-mail: Gregg . Wade@rmc . ca 
* Based on observations obtained using the MuSiCoS spectropo- 
larimeter at the Pic du Midi observatory, France 



namo that has since turned off . As there is no particular reason 
to suspect that similar processes should not occur during the 
formation of more massive stars, it seems reasonable a priori 
to expect that magnetic fields of similar structure and intensity 
should also exist in O-type stars. 

In the absence of direct magnetic measurements, this view 
has been supported by the wide-spread occurrence of variabil- 
ity in the winds of O-type stars, which may provide indirect 
evidence for the presence of dynamically important magnetic 

1 Credible quantitative models invoking contemporaneous dy- 
namos operating in the convective core have also been proposed (e.g. 
Charbonneau & MacGregor 2001), but these models generally have 
significant difficulty explaining the intensities and topologies of the 
observed fields, their diversity, lack of correlation of field with angu- 
lar velocity, as well as the young ages of many magnetic stars. 
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fields in their atmospheres. During the past 20 years, sustained 
optical and UV spectroscopic observations have shown that the 
winds of O-type stars are highly structured, exhibiting both 
coherent and stochastic behaviour, as well as cyclical variabil- 
ity on timescales ranging from hours to days. 2 A key result of 
this work has been the conclusion that a major component of 
this variability results from rotational modulation of structures 
imposed on the wind by some deep-seated process. Magnetic 
fields are pres ently considered a likely source of such structures 
(e.g. JCranmer & Owockifl99(J) 3 . 

Although by no means prototypical, 9 l Orionis C 
(HD 37022; HR 1895) is perhaps the best example of an 
O-type star with distinctive, periodic variability of its spec- 
troscopic stellar wind features. It is a very young, pecu- 
liar ~07 star, and the brightest and hottest member of the 
Orion Nebula Cluster. It exhibits strictly periodic spectroscopic 
variability which is strongly suggestive of a magnetic rota- 
tor: Ha and Hen /14686 emission, peculiar ultraviolet Civ 
/LU548, 1550 wind lines, photospheric absorption lines, and 
ROSAT X-ray emission all appear to vary with a single well- 
defined period of 15-422 + 0.002 d {stahl et alJll993l Il996t 
IWalborn & Nichols! 1994ldagne et alll997t) . The similarity of 
this beha viour to that of some magnetic A- and B-type stars 
(see, e.g., Ish ore & Brown 1990) has motivated the suggestion 
first made bv lStahl et alJ i 19961) that 6 1 Ori C also hosts a fos- 
sil magnetic field which confines the stellar wind and produces 
the observed variability by rotational modulation. Such a phe- 
nomenon appears to have been first suggested (for the O star 
£ Pup) by Moffat & Michaud (1981) (although no field has ever 
been detected in that star). 

In order to test this suggestion. lDonati & Wadel i 19991) ob- 
tained longitudinal magnetic field measurements of 9 1 Ori C, 
but failed to detect the prese nce of a photosph eric field. 
However, in a seminal paper, iDonati et alJ (J2002) reported 
the detection of a fossil magnetic field in the photosphere of 
6 l Ori C based on 5 measurements of the Stokes V profiles 
of selected photospheric absorption lines. This detection pro- 
vided the first empirical support that dynamically important 
magnetic fields exist in O-type stars, and that their behaviour 
is directly linked to spectroscopic variability. 4 It also repre- 
sents the youngest main sequence star in which a fossil-type 
field has been detected (comparable in age to NGC 2244-334; 
Bag nulo et all2004l) . Since 6 1 Ori C is clearly a pivotal object, 
it is important that this detection be corroborated and the deter- 
mination of its magnetic properties be refined. 

The primary goal of the present paper is to report confir- 
mation of the detection of a magnetic field in the photospheric 
lines of 9 l Ori C. Our new spectropolarimetric observations, 

2 See,e.g.,Fullertonl2003) for a recent review of the characteristics 
of variability in hot-stars winds. 

3 Magnetic fields are also seen by some investigators as an ingredi- 
ent necessary to explain intrinsic X-ray fluxes and non-thermal radio 
emission fro m massive star s. 

4 Recently, Don ati et all 12005) reported the detection of a 1.5 kG 
dipolar magnetic field in the O-type spectrum variable HD 191612 
[Of?p], which exhibits stellar wind variations with a period of 538 
days. They suggested that HD 191612 represents an evolved version 
of 6 1 OriC. 



which provide a substantially larger data set that samples the 
15 < ?422 rotational period more fully, are described in §2. The 
analysis of the Least-Squares Deconvolved circular polarisa- 
tion spectra, both in terms of the mean longitudinal magnetic 
field, (B z ), and the LSD mean Stokes V profiles, is described 
in §3. In §4 we discuss the variability characteristics of fea- 
tures in the Stokes / spectrum, and comment on consistency 
with both the mag netically-confined wind shock model of 
IDonati et alJ d2002t) as well as ver y recent MHD simulat i ons of 
magnetically-channelled winds by ud-Doula & Owocki (2002) 
and iGagne et alJ J2005I) . Finally, in §5 we summarise our re- 
sults, and discuss implications for our understanding of the ori- 
gin and evolution of magnetic fields in intermediate- and high- 
mass stars, and for our understanding of wind modulation in 
massive stars. 

2. Observations 

Circular polarisation (Stokes V) spectra of 8 1 Ori C were ob- 
tained during the period 1997-2000 using the MuSiCoS spec- 
tropolarimeter mounted on the 2 metre Bernard Lyot tele- 
scope at Pic du Midi observatory. The spectropolarimeter 
consists of a dedicated polarimetric module mounted at the 
Cassegrain focus of the telescope and connected by a dou- 
ble optical fibre (one for each orthogonal polarisation state) 
to the table-mounted cross-dispersed echelle spectrograph. The 
spectrogr aph and pola rimeter module are described in d etail 
bv lBaudrand & Bohml i 19921) and bv lDonati etaD dl999l) . re- 
spectively. The standard instrumental configuration allows for 
the acquisition of circular or linear polarisation spectra with a 
resolving power of about 35,000 throughout the range 4500- 
6600 A. 

A complete circular polarisation observation consists of 
a series of 4 sub-exposures between which the polarimeter 
quarter-wave plate is rotated back and forth between position 
angles (of the plate fast axis with respect to the polarising 
beamsplitter fast axis) of -45° and +45°. This procedure results 
in exchanging the orthogonally polarised beams throughout the 
entire instrument, which makes it possible to reduce system- 
atic errors (due to interference and other effects in the tele- 
scope and polarisation optics, instrumental drifts, astrophysical 
variability, etc.) in spectral line polarisation measurements of 
sharp -lined stars to below a level of about 10~ 4 dWade et alJ 
2000). 

In total, 45 Stokes V spectra of 1 Ori C were obtained over 
4 observing runs in 1997 February, 1998 February/March, 2000 
February/March, and 2000/2001 December/January, with peak 
signal-to-noise ratios (S/N) of typically 250 per pixel in the 
continuum. The 6 spectra obtained during 1997 and 1998 have 
already been discussed by D onati & Wadel \ 19991) . Some fring- 
ing is visible in the Stokes / spectra (maximum amplitude ~ 
1 % peak-to-peak in the red), although no fringing is evident in 
Stokes V. Such effects are present in most polarimeters, and 
generally result from internal reflections producing secondary 
beams coherent with the incident beam, but with important 
(wavelength dependent) phase differences (Semel 2003). This 
fringing limits the effective S/N of the Stokes / spectra to about 
150:1 at Ho-, and to about 350:1 at He n ,14686. 
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Table 1. Journal of obs ervations. Phases a re calculated accord- 
ing to the ephemeris of IStahl et al] Jl9961) . JD = 2448833.0 + 
15.422 E. The peak S/N per pixel in the continuum is quoted in 
the final column. 



Table 2. Binned LSD Stokes V profiles and longitudinal mag 
netic field measurements 



No. 


Date 


HJD 


Phase 


t 

' i_'Xp 


S/N 






(-2450000) 




[minutes] 


pixel ^ 


01 


1997 Feb 20 


0500.3453 


0.1147 


40 


280 


02 


1997 Feb 22 


0502.3429 


0.2442 


40 


450 


03 


1997 Feb 23 


0503.3230 


0.3078 


40 


140 


04 


1997 Feb 25 


0505.3160 


0.4370 


40 


300 


Uj 


1 QQC CoK 1 ^ 


UOOU.4UZ1 


n a ^1 a 
U.4010 


-in 
4U 


o^n 

zjU 


uo 


1 QQC CoK OA 

iVVo reD ZD 


Uo / 1. jjj4 


U. 1 /UO 


-in 
4U 


07n 

Z /U 


07 


2000 Feb 03 


1578.4219 


0.0198 


40 


250 


08 


2000 Feb 03 


1578.4536 


0.0218 


40 


240 


09 


2000 Feb 04 


1579.3932 


0.0828 


40 


250 


10 


2000 Feb 04 


1579.4235 


0.0847 


40 


270 


11 


2000 Feb 04 


1579.4539 


0.0867 


40 


320 


12 


2000 Feb 12 


1587.3739 


0.6003 


40 


340 


13 


2000 Feb 15 


1590.3822 


0.7953 


40 


310 


14 


2000 Feb 15 


1590.4131 


0.7973 


40 


260 


15 


2000 Feb 21 


1596.3907 


0.1849 


40 


270 


16 


2000 Feb 21 


1596.4211 


0.1869 


40 


370 


17 


2000 Feb 21 


1596.4518 


0.1889 


40 


350 


18 


2000 Feb 22 


1597.3515 


0.2472 


40 


410 


19 


2000 Feb 22 


1597.3823 


0.2492 


40 


400 


20 


2000 Feb 22 


1597.4133 


0.2512 


40 


410 


21 


2000 Feb 22 


1597.4439 


0.2532 


40 


380 


22 


2000 Feb 24 


1599.3534 


0.3770 


40 


410 


23 


2000 Feb 24 


1599.3840 


0.3790 


40 


400 


24 


2000 Feb 24 


1599.4156 


0.3811 


40 


340 


25 


2000 Feb 25 


1600.3465 


0.4414 


40 


410 


26 


2000 Feb 25 


1600.3780 


0.4435 


40 


410 


27 


2000 Feb 25 


1600.4105 


0.4456 


40 


390 


28 


2000 Feb 26 


1601.3610 


0.5072 


40 


340 


29 


2000 Feb 26 


1601.3922 


0.5092 


40 


340 


30 


2000 Feb 26 


1601.4229 


0.5112 


40 


310 


31 


2000 Feb 27 


1602.3490 


0.5713 


40 


360 


32 


2000 Feb 27 


1602.3881 


0.5738 


40 


350 


33 


2000 Feb 27 


1602.4191 


0.5758 


40 


310 


34 


2000 Feb 28 


1603.3767 


0.6379 


40 


240 


35 


2000 Feb 28 


1603.4077 


0.6399 


40 


300 


36 


2000 Feb 28 


1603.4391 


0.6420 


40 


230 


37 


2000 Mar 02 


1606.3531 


0.8309 


40 


340 


38 


2000 Mar 02 


1606.3844 


0.8329 


40 


390 


39 


2000 Mar 02 


1606.4161 


0.8350 


40 


340 


40 


2000 Mar 04 


1608.3696 


0.9617 


40 


470 


41 


2000 Mar 04 


1608.4005 


0.9637 


40 


410 


42 


2000 Mar 04 


1608.4311 


0.9657 


40 


360 


43 


2000 Mar 05 


1609.3605 


0.0259 


40 


350 


44 


2000 Mar 05 


1609.3910 


0.0279 


40 


310 


45 


2000 Mar 05 


1609.4216 


0.0299 


40 


260 



Mean 


Binned 


(B : ) ± erg 


N LSD 




Phase 


Spectra 


(G) 


(%) 




0.0550 


01,07-09,10,11,43-45 


480 ± 95 


3.97 10" 


-2 


0.1829 


06,15-17 


590 ± 133 


5.38 10" 


-2 


0.2491 


02,18-21 


334 ± 96 


3.74 10" 


•2 


0.3613 


03,22-24 


120 ± 130 


5.04 10" 


-2 


0.4459 


04,05,25-27 


-110 ± 118 


4.27 10" 


-2 


0.5415 


28-33 


130 ± 110 


4.1210" 


-4 


0.6301 


12,34-36 


-96 ± 164 


5.99 10" 


-2 


0.7964 


13,14 


14 ± 226 


8.41 10" 


-2 


0.8330 


37-39 


340 ± 153 


5.67 10' 


•2 


0.9637 


40-42 


589 ± 116 


4.89 10 


•2 



In addition to exposures of 6 l Ori C, during each run obser- 
vations of various magnetic and non-magnetic standard stars 
were obtained which confir m the nomi nal operation of the in- 
strument (see, e.g., IShorlin et al.ll200 2). The observing log is 
shown in Table 1 . 



3. Least-Squares Deconvolution 

As demonstrated by, e.g., iDonati & Wadel Jl999), the Least- 
Squares Deconvolution (LSD) multi-line procedure can pro- 
vide enormous improvement in the precision of spectral line 
polarisation measurements as compared with individual line 
measurements. We be gan using the 13 lines employed by 
IDonati & Wadel dl999l) for the LSD analysis. After exploring 
the effect of removing lines from this list, it became clear based 
on the shape of the Stokes / profile that only 3 of the lines in 
the mask were contributing relatively uncontaminated photo- 
spheric profiles. In the end, an LSD mask including only Oni 
/i5592, Civ /15801 and Civ /15811 was used, with a mean wave- 
length A = 5730 Aand a mean Lande factor z = 1.18. We 
found that the resultant LSD S/N was dependent on the ve- 
locity bin size of the extracted LSD profiles. Ultimately, ex- 
periment yielded a best S/N for an bin size of 13.5 kms" 1 . 
Therefore, each pixel in the extracted Stokes / and V LSD pro- 
files corresponds to a velocity interval of 13.5 kms -1 . Each 
extracted LSD profile was determined from 4 individual pho- 
tospheric features (Civ /15801 appears in two separate orders 
in each spectrum). 

To further reduce the noise, LSD profi les were binned in 
phase, according to the 15.422 day period o flStahl et alJ (fl996). 
weighting each pixel according to the inverse of its squared 
error bar. The final binned profiles, which will be used for all 
further LSD analysis, are summarised in Table 2. 

4. Magnetic field geometry 

Following the approach of Dona ti"et"al] d200ll) . the best-fit ge- 
ometry of the dipolar magnetic field was determined using two 
complementary procedures. 

First we modeled the variation of the mean longitudi- 
nal magnetic field. The longitudinal field (B z ) and its asso- 
ciated error cr B were inferr ed from each set of LSD profiles, 
in the manner describ ed bv IDonati et alJ d 19971) and modified 
bv IWade et al. (2000). Measurements were made in the range 
+80 kms" 1 around line centre, in both the Stokes V and diag- 
nostic N LSD profiles. The N, or null profiles are calculated 
from analysis of Stokes V CCD frames obtained at identical 



4 



G.A. Wade et al.: The magnetic field and confined wind of the O star 9 [ Orionis C 




_400 ' 1 1 1 1 1 1 1 1 1 l 

0.5 1 

PHASE (P= 15.422 d) 



Fig. 1. Longitudinal magnetic field variation (in G) of 6 1 Ori C 
phased according to the ephemeris described in the text. Upper 
panel - measurements from Stokes V profiles; Lower panel - 
measured from diagnostic N (null) profiles 

waveplat e angles, and a re nominally consistent with zero (see, 
e.g., Don ati et all 1 997). The N profiles provide a powerful di- 
agnosis of systematic errors in the polarisation data. The mea- 
surements are summarised in Table 2 and are shown, phased ac- 
cording to the 15.422 d period, in Fig.^ The reduced^ 2 of the 
(phased) Stokes V measurements of (B z ) is 1.0 for a first-order 
sine fit (fitting zero-point, amplitude and phase), and 10.2 for 
the null-field hypothesis (i.e. a straight line through (B.) = 0. 
The reduced^ 2 of the N (B z ) measurements is 0.8 for a first- 
order sine fit (fitting zero-point, amplitude and phase), and 1 .2 
for the null-field hypothesis. Based on these results, for the 10 
longitudinal field data discussed here, the null-field hypothe- 
sis can be confidently ruled out for the Stokes V (B : ) variation 
(false-alarm probability <K 0. 1 %), whereas a null magnetic field 
is consistent with the N measurements within 2cr. 

The l st -order sine fit to the Stokes V (B z ) data is charac- 
terised by a maximum of 550 + 85 G, a minimum of -35 + 85 G 
and a phase of maximum cpa = 0.05 + 0.05 (all uncertainties 
lcr). 

Then, using a modified version of the programme Fldcurv, 
we calculated model longitudinal field variations correspond- 



ing to a large number of dipolar magnetic field configurations, 
varying the dipole intensity B d as well as the obliquity angle 

B. Fldcurv computes the surface distribution of magnetic field 
corresponding to a particular magnetic geometry, then weights 
and integrates the longitudinal component over the hemisphere 
of the star visible at a given phase (for this procedure we have 
assumed a limb-darkening coefficient of 0.4). 

In order to uniquely identify the magnetic geometry, we 
must specify four parameters: The magnetic dipole polar 
strength B d , the inclination of the stellar rotational axis to the 
observer's line-of-sight i, the obliquity of the magnetic axis 
with respect to the rotational axis /3, and the phase of closest 
passage of the positive magnetic pole to the line-of-sight, <f>o. 
The phase variation of the longitudinal field allows us to con- 
strain 3 of these parameters; typically B d ,B and 0o are deter- 
mined by fitting the magnetic c urve, w hereas i is constrained 
using other data, [Donati et al. (2002) derive i = 45 ± 20° 
based on the rotational period and inferred vsim of G l Ori 

C, as well as the characteristics of the optical and UV spec- 
troscopic variability. Their conclusions are independ ently sup- 
ported by the optical spectroscopic study of ISimon-Diaz et all 
(2005), who find i = 44 ± 12°, and qualitatively by the orig- 
inal Magnetically-Confined Wind Shock modelli ng of th e X- 
ray variation of this star by iBahel A Montr^eTkl dl 997hh . We 
therefore adopt i — 45 + 20° for the modelling performed in 
this paper. 5 

For each computed model (over 40,000 in total), the re- 
duced x 1 of the model fit to the observations was calculated. 
The 'map' of these ^ 2 s, shown in the upper frame of Fig. |2 
(where black regions represent acceptable models, whereas 
white regions represent models rejectable at more than the 95% 
[approximately 2<x] confidence level 6 ), indicates that for a ro- 
tational axis inclination i - 45°, acceptable magnetic models 
are characterised by B d = 1900 + 250 G and B = 50 + 8° 
(where associated uncertainties are quoted at lcr). If we allow 
a 20° uncertainty in our assumed value of i, we obtain Bd = 
2400 + 300 G and B = 30 + 8° (z = 65°) and B d = 2600 + 450 G 
and B — 67 + 5° (;' = 25°). Therefore our modeling of the longi- 
tudinal field variation constrains the dipole magnetic field ge- 
ometry of l Ori C to 1650 <B d < 3050 G and 22° < B < 73°. 

Our second procedure involves direct modeling of the LSD 
Stokes / and V profiles, in a manner similar to that described 
first by Donat i et alJ ( 1200 II) . In particular, all line profile mod- 
eling was accomplished using the LTE polarised synthesis code 
Zeeman2 (Landstreet 1988; Wade et al. 2001) and assuming an 
Atlas9 model atmosphere with T e ff = 40, 000 K and log g 
4.0. 

We began by finding a phase-independent model fit to the 
Stokes / profiles, varying the line equivalent width, projected 

5 Although all available evidence points to the adopted inclination, 
we have performed additional modelling to explore the sensitivity of 
more extreme inclination angles on the derived magnetic field geom- 
etry. For inclinations as small as 10° and as large as 80°, we find that 
the conclusions of this paper are not qualitatively affected. 

6 Confidence interv als were obtained using x 1 probability tables 
from lBevingtonl <1969t) . According to these tables, the 2cr (95%) confi- 
dence interval, considering each model parameter independently, cor- 
responds to an increase in totals 2 of 3.84. 
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rotational velocity v sin i, and microturbulent broadenin g C as 
free parameters. Assuming vsin/ = 20 kms 1 JPonati et alJ 
2002), we found a microturbulent broadening £ = 30 kms -1 
(corresponding to a total thermal +turbulent broadening of 
abou t 3 1 km s~' ) provi ded a best-fit to the profiles. As reported 
bv lDonatietai]j2002l) . the LSD Stokes / profiles vary system- 
atically according to rotational phase, with a maximum depth 
of about 12% of the continuum at phase 0.0, and a minimum 
depth of about 9% of the continuum at phase 0.5. This varia- 
tion is quite probably related to variable contamination of the 
photospheric profiles by the wind spectrum; see §5. 

The next step was to calculate Stokes V profiles corre- 
sponding to a large number of dipole surface magnetic field 
configurations, varying the dipole parameters B and B d for se- 
lected values of i. Finally, we compared each of these calcu- 
lated profiles (about 28,000 in all) with the observed Stokes V 
profile, calculating the reduced* 2 and building up a "map" of 
their agreement as well. 

The results are illustrated in the lower frame of Fig. [2] (for 
i = 45°) and in Fig.|3(for i = 25° and i = 65°) for 95% (2cr) 
confidence. For i = 45°, we obtain a somewhat smaller value 
for the intensity of the best-fit dipole, B d = 1300 ± 150 G, 
and B = 50° ± 6° (again, all quoted uncertainties are 1(f). For 
i = 25° we obtain B d = 1350 ± 150 G and B = 62° ± 6°, while 
for i = 65° we obtain B d = 1600 ± 200 G and B - 33° + 6°. 
Therefore our modeling of the LSD mean Stokes profile varia- 
tion constrains the dipole magnetic field geometry of 9 l Ori C 
to 1 150 < B d < 1800 G and 27° < B £ 68°. The reduced* 2 for 
the null-field hypothesis for the V measurements is 1.11, ver- 
sus 0.76 for the best-fit dipole model, indicating that the null- 
field hypothesis can be confidently ruled out for the Stokes V 
spectrum (false-alarm probability <K 0.1%). At the same time, a 
similar analysis of the N diagnostic null LSD profiles provides 
a best-fit model of B d = 150 ± 225 G with reduced x 2 of 0.69, 
consistent with a null field and indicating that no significant 
fringing contamination of the N profiles exists. 

The low reduced y 2 statistics obtained for both Stokes V 
and N suggest that the error bars are slightly overestimate d 
(by about 10%). Based on the r esults of Wa de et alJ J2000), 
this is not unexpected. Moreover, the difference in best-fit re- 
duced x 2 obtained for the V versus the N profiles (0.76 ver- 
sus 0.69) is marginally significant, indicating that (not surpris- 
ingly) the Stokes V signatures are not quite fit to within their 
error bars, and suggesting that other unmodelled effects may 
be at play. The computed Stokes profiles corresponding to the 
best-fit model are compared with the mean LSD profiles in Fig. 
3. 

The general solutions for the magnetic field geometry ob- 
tained using these two methods are in good agreement. The B d 
solutions for i = 45° differ at the 1.5cr level. Although not es- 
pecially significant, this difference could be attributable to the 
more approximate weighting of the local field in the Fldcurv 
model, or possibly contributions to the amplitude and variabil- 
ity of the Stokes V signatures by emission contamination of the 
photospheric profiles by the circumstellar material. 

Due to the more sophisticated nature of the direct Stokes 
V fitting procedure, we adopt these results as our formal con- 
straints on the surface magnetic field of 9 1 Ori C. These re- 



sults a re perfectly consistent with those derived bv lPonati et alJ 
il2002l) . but are to be preferred since they are derived from a 
larger sample of spectra. We therefore confirm the detection 
and the ge ometrical cha r acteris tics of the magnetic field dis- 
covered by Donati et al] J2002I) . Moreover, we can affirm the 
sinusoidal nature of the longitudinal field variation and con- 
fidently define the phases of maximum and minimum (B z ) 
(0.05 ± 0.05 and 0.55 ± 0-05, respectively) according to the 
ephemeris o f Stah l et alJJl996l) . 

The geometry of the magnetic field an d circumstellar ma- 
terial of l Ori C is sketche d in Fig. 1 of ISmifh & Fullertoi] 
(2005) or Gag ne et alJ J2005t) . and is illustrated as a series of 
animations at www. astro .udel . edu/tloc. 

5. Line-profile variations 

Phase-resolved variations of the Stokes / profiles from the 2000 
time series are displayed in an image format for selected transi- 
tions in Figures|5HEl in these images, which are commonly re- 
ferred to as "dynamic spectra", individual spectra occupy hor- 
izontal strips, which are stacked vertically to show systematic 
variations with time. The transitions are grouped to illustrate 
variations that are predominantly circumstellar (Fig. [5] and [6j 
and photospheric (Fig. 0, with selected He 1 lines serving as 
hybrid cases (Fig.[S}. 

The line-profile variations of 6 l Ori C have been previously 
illustrated as dynamic spectra by ISlahletalJ ll99fl for Ha, 
Hen ,14686, C iv AA1548, 1550, Si rv AA 1393, 1402, and Ora 
,15592) and bv lReiners et ail J2000l for He 1 A4 47 1 and ,14713, 
Civ A581L and O m A5592v lstahl et alJJl996h represented the 
variations as differences with respect to either (a) line profiles 
from a similar star [15 Monocerotis; spectral type 07 V((f))] 
that are at most modestly contaminated by stellar-wind emis- 
sion; or (b) minimum absorption profiles from the time se- 
ries. Reiners et al. (2000) illustrated the variations in the spec- 
tra themselves, i.e., without renormalizing by a template spec- 
trum. In contrast, the dynamic spectra illustrated here represent 
variations with respect to the mean spectrum, which is reason- 
ably well defined by the approximately uniform sampling of 
this time series. The advantage of displaying quotient spectra 
with respect to the mean profile is that the relative amplitude 
of variations can be compared in a meaningful way between 
(unsaturated) lines of different shape and strength. Whatever 
template is used to normalize a time series, it is worth remem- 
bering that the emergent spectrum is a very nonlinear function 
of, e.g., hydrodynamic variables such as density. Consequently, 
quantitative interpretation of the line profile variations in terms 
of physical parameters generally requires comparison with a 
model. 

5.1. Circumstellar lines 

Figures|5]and|S]illustrate the line-profile variations for Ha, HB, 
and He n ,14686, all of which are dominated by stellar wind ma- 
terial in the magnetosphere. All three lines exhibit qualitatively 
similar patterns of variability. In Fig. [5] the contrast is adjusted 
separately for each line in order that the entire variation can be 
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seen, while in Fig. [5] a fixed high-contrast "stretch" is used to 
enhance the visibility of weaker variations. 

Fig-IUconfirms that the basic modulation consists of a broad 
emission excess that attains maximum strength at phase 0.0 
(when the magnetic equator is viewed approximately face-on), 
decreases to a minimum near phase 0.5 (when the magnetic 
equator is viewed approximately edge-on), and increases again 
starting near phase 0.75. Although the overall symmetry of the 
emission feature is difficult to determine due to contamination 
from the nebular component in Ha and H/3, it is blue-shifted 
near phase by perhaps as much as 100 kms~'. As the emis- 
sion fades, the primary variations are more nearly centered on 
the rest velocity of the star. 

However, Fig. [6] also indicates the presence of a second 
component, which is particularly prominent in Hen /J4686. It 
consists of an approximately stationary excess with respect to 
the template at v r * 200 kms~' between phases ~0.4 and 0.6, 
which achieves maximum strength near phase 0.5 (i.e., when 
the magnetic equator is viewed edge-on). Although it is often 
difficult to determine whether an increase in the quotient spec- 
trum represents excess emission or reduced absorption with re- 
spect to the template spectrum, it is clear from examination 
of the line profiles themselves that this component is due to 
emis sion. This com ponent has been noted previously in Ha 
by IStahl et al.l ll 19931) . who concluded that it was responsible 
for the smaller minimum near phase 0.5 in the "M-shaped " 
equivalent width curve; see, e.g., Fig. 4 of IStahl et alJ Jl996l) . 
However, the origin of this feature was not discussed. 



5.2. Photospheric lines 

Fig.[7]shows the line-profile variations of "photospheric" lines, 
i.e., lines that are not obviously contaminated by emission from 
the stellar wind or material trapped in the magnetosphere. The 
lines included in this category are Civ /L15801, 5811, Om 
/15592, and He n A541 1 . The first three of these are combined to 
obtain the LSD measurements of the magnetic field of 6 l Ori C. 

The variations in all these lines are qualitatively similar, in- 
cluding Civ ,45801, which is not illustrated. They are deeper 
with respect to the mean profile between phases 0.0 and ~0. 15; 
less deep between phases 0.2 and 0.6 (possibly later); and 
deep er again betwee n phases ~0.8 and 1 .0. As previously noted 
bv lStahl et alJ Jl996l) . the variations appear to move from blue- 
to-red (i.e., in the sense of rotation), particularly near phase 0, 
but are confined to the central region of the line profile. The be- 
haviour illustrated here is very similar to that exhibited in other 
dyna mic spectra, despi te di fferences in the temp lates used; see, 
e.g. JStahl et aDdl996h and lReiners e7aTIJ2000l) . 

Although not illustrated here, several other high-excitation 
lines were also investigated. No significant variations were vis- 
ible in the Hen /14542 and Siiv /14654 absorption lines or 
the Cm /I5696 selective emission feature. A weak modulation 
might be present in Si iv A463 1 . 



5.3. He i lines 

The line profile variations of the He i lines illustrated in Fig. [8] 
span the behaviour of "wind" and "photospheric" lines. In par- 
ticular, the strong Hei /15876 triplet exhibits "circumstellar" 
variations. This is not surprising, since the sensitivity of Hei 
/15876 to low-density environments is a well-known conse- 
quence of non-LTE physics. However, the "circumstellar" vari- 
ations in Hei /15876 occur at smaller velocities than indicated 
in Fig. |3 

In contrast to He i /15876, the weaker He i /14713 triplet and 
/14921 singlet show "photospheric" variations that are very sim- 
ilar in phasing and character to those illustrated in Fig. [7] The 
variations in Hei ^47 13 are substantially weaker than those 
in Hei /14921, even though the lines themselves have similar 
strength. Although not illustrated here, the Hei /15015 singlet 
also shows weak "photospheric" variations. 

5.4. Origin of the line profile variations 

The primary modulatio n of the circumstellar line profiles has 
been modelled by Do nati et"al] J2002I) in terms of the "mag- 
netically confined wind shock model" of iBabel & Montmerld 
( 1997b a). The key feature of this model is that a sufficiently 
strong, dipolar magnetic field channels wind material to the 
magnetic equator, where it collides with material similarly 
channeled from the opposite hemisphere. The consequences in- 
clude substantial shock heating of the gas, the development of 
extended cooling zones, and the creation of a zone of enhanced 
density at the magnetic equator. With this model and their re- 
vised mappin g between ro t ationa l and magnetic phase (which 
we confirm), Donati et 

ID £2002) were able to reproduce the 
fundamental features of the primary modulation; see, e.g., their 
Fig. 12. As a result, the maximum emission from circumstel- 
lar material is now understood to come from the cooling region 
above and below the magnetic equator, which is viewed face- 
on at phase 0.0. In contrast, the observer views the magnetic 
equatorial plane edge-on at phase 0.5, at which time the cool- 
ing disk is viewed edge-on against the stellar photosphere and 
its contributions t o the Ha prof i le are minimized. Despite the 
generally good fit. iDonati et alJ J2002) noted some difficulties 
with their model profiles, which also did not reproduce the red- 
shifted emissi on feature near pha se 0.5 

Recently, Gag ne etall |2P05) presented a magnetohydro- 
dynamic (MHD) model of the wind of 9 l Ori C that confirms 
and extends the basic picture provided by the "magnetically 
confined wind shock" model. The new "magnetically chan- 
neled wind shock" model follows the same com putational ap- 
proach described bv bd-Doula & Owo cki (2002) by allowing 
for the dynamic competition between the outward forces of 
radiation pressure and the channeling and confinement of the 
magnetic field. However, in addition to allowing for the feed- 
back of the outflow on the geometry of the magnetic field, the 
newer model also incorporates a detailed treatment of the en- 
ergy balance in the wind, including the effects of compressive 
heating d ue to the strong sho cks in the vicinity of the magnetic 
equator. Gag ne et alJ J2005) show that this model successfully 
explains both the periodic modulation of X-ray emission exhib- 
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ited by 9 1 Ori C and the basic features of the X-ray emission 
lines. 

The new MHD model also suggests that the circumstellar 
environment of 6 1 Ori C is much more dynamic than previ- 
ously supposed. In particular, the models show that the amount 
of material confined to the region of the cooling disk is limited 
by outflow through the disk beyond the Alfven radius (i.e., the 
point at which the magnetic field can contain material chan- 
neled to the region of the magnetic equator) as well as in- 
fall in the inner regions. This infall occurs sporadically when 
cool, compressed material in the magnetic equator becomes too 
dense to be supported by radiative driving, and consequently 
falls back onto the stellar surface along distorted field lines. 
Althoug h much o f the infalling material is too cool to emit X- 
ravs. lGagne et alJ (2005) found evidence that the radial veloci- 
ties of X-ray emission lines were systematically red-shifted by 
as much as 93 kms -1 when the disk was viewed edge on (i.e., 
ne ar phase 0.5) 

ISmith & FuUertorJll2005l) also recognized that the presence 
of infalling material near the magnetic equator provides a clue 
to the origin of the red-shifted emission feature near phase 0.5 
in line profiles dominated by circumstellar material (Fig. [6}. 
They attributed this feature to a column of optically thick ma- 
terial flowing inward at the magnetic equator, and used this 
interpretation to explain previously unnoticed modulations of 
the Civ and Nv resonance doublets at modest red-shifts of 
250 kms~' or less. We similarly attribute the red-shifted emis- 
sion features near phase 0.5 in Fig.[6](and also for He i /15876 
in Fig. |HJ to the infalling material seen in the MHD models. In 
this context, the fact that the emission feature is seen at smaller 
red-shifted velocities in He i /15876 (~50 km s ; possibly over 
a smaller range of phases) than its counterpart in He n /14686 
(~200 kms -1 ) suggests that He + is formed closer to the radius 
where infall is initiated, and that the material heats up as it ac- 
celerates on its return to the star. 

Although the presence of infalling material appears to be 
a robust consequence of magnetic channeling for a star like 
6 l Ori C, the current generation of MHD simulations indicates 
that it is an occasional occurrence, which is not tied to any 
specific rotational or magnetic phase. Consequently, the con- 
sistent presence of the red-shifted features near phase 0.5 in 
time series of optical, UV, and X-ray lines obtained over many 
rotational cycles is surprising. It remains to be seen whether re- 
fined MHD models can reproduce the approximate steady-state 
of infall implied by these observations. Since the mass-loss rate 
of 6 l Ori C is poorly constrained, one straightforward (though 
arbitrary) solution might be to increase the value assumed in 
the MHD simulations. An increase in the basal mass flux may 
slightly alter the magnetic geometry above the stellar surface, 
but it will certainly increase the rate at which material stag- 
nates near the magnetic equator, perhaps to the point where the 
frequency of infall in the simulations matches the high duty- 
cycle implied by the observations. However, a similar steady- 
state accretion also seems to be required for the magnetic B star 
/3 Cep, for which the mass-loss rate is well constrained. Such 
a scenario therefore does not appear to be applicable to f3 Cep, 
but it is not obvious that this conclusion can be extended to the 
case ofO 1 OriC. 



The red circumstellar emission component at phase 0.5 
may also p l ay a r ole in the variations of "p hotospheric" lines. 
IStahl et alJ ill 99(h and iReiners et all i2QQ(ii interpreted these 
variations in terms of an absorption excess at phase 0.0, be- 
cause the alternative - an emission excess at phase 0.5 - con- 
tradicted the behaviour of the primary modulation observed 
in emission l ines, which achieves maximum at phase 0.0. 
Rei ners et alJ J2000h modelled the photospheric variations in 
terms of asymmetrical distributions of spots characterized by 
(a) reduced abundance near the magnetic poles and (b) en- 
hanced abundance around the magnetic equator. They found 
reasonable agreement in both cases, but particularly for (b). 
Unfortunately, the mappi ng between rotational p hase and mag- 
netic g eometry used by IStahl et aD i 19961) and IReiners et alJ 
(2000) turned out to be incorrect (see iDonati et alJl2002l) . so 
their interpretations require revision. 

More recently, ISimon-Diaz et alJ J2005I) attempted to ex- 
plain the variations of the photospheric lines as a consequence 
of additional continuum light from the disk-like structure at 
the magnetic equator. However, they also assumed that rota- 
tional phase 0.0 corresponds to the configuration when the 
magnetic equator is viewed edge-on, in order that the maxi- 
mum contamination (hence minimum line depth and equiva- 
lent width) will be seen at rotational phase 0.5. This mapping 
between rotational phase and magnetic geometry is not sup- 
ported by the magnetic field measurements. Furthermore, con- 
tinuum variations on the rotational period of the size they pre- 
dict (0.16 magnitu des, peak-to-peak) a re excluded by the pho- 
tometry of Ivan Genderen et al. I (U985J. Consequently, the ex- 
pl anation for the photosphe ric line-profile variations proposed 
by Simo n^Diaz et alJ J2005I) is not viable. 

Instead, we note that the photospheric line profiles are less 
deep at the same time that the circumstellar lines exhibit the 
red emission component. This behaviour is seen particularly 
well in Fig. [8] which shows that the red-wing emission of 
the "circumstellar" Hei /15876 line is in phase with the mid- 
cycle absorption minima in the other, "photospheric" He i lines. 
Consequently, if the explanation for the red-shifted component 
at phase 0.5 in circumstellar lines in terms of infall is correct, 
then the synchronized variations of the "photospheric" lines 
might also be attributable to the infalling material. From this 
perspective, the change in photospheric line depth is due to par- 
tial emission filling near phase 0.5, rather than excess absorp- 
tion. Quantitative modelling based on magnetohydrodynamic 
simulations is required to determine whether the geometry and 
properties of the infalling material can affect high-excitation 
lines like Civ /L15801, 5811, e.g., by producing an optically 
thick screen that shadows a small fraction of the stellar disk 
(thereby minimizing continuum photometric variations) but is 
sufficient to change the appearance of the red half of a photo- 
spheric line profile. In particular, it remains to be seen whether 
sufficient density accumulates near the stagnation point to ex- 
plain the presence of high-excitation emission at small infall 
velocities. 
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6. Summary and conclusions 

Using a new series of 45 Stokes / and V spectra obtained 
with the MuSiCoS spectropolarimeter at Pic du Midi obser- 
vatory, we have detected the photospheric magnetic field of 
the young 07 Trapezium memb er 6 1 Ori C. We confirm and 
extend the conclusions of Donati et"afl J2002I) : that the Stokes 
V variations are consistent with a dipolar magnetic field with 
a polar strength between 1150 and 1800 G, and an obliquity 
27° < /3 < 68° if i = 45 + 20°. Moreover, we demonstrate 
that the variation of the longitudinal magnetic field is sinu- 
soidal to within the errors; the phase of maximum longitudinal 
field i s 0.05 + 0.05 according to the ephemeris of IStahl et alJ 
and the longitudinal field varies between —35 ± 85 G 
and +550 + 85 G. We have also exploited our high-resolution 
Stokes / spectra to study the cyclical variations of spectral 
absorption and emission lines formed in the photosphere and 
wind of 9 l Ori C. We confirm th e varia bility properties reported 
previously by, e.g.. lStahl et alJ {l<)96), and highlight evidence 
that suggests the presence of infalling material in the magnetic 
equatorial plane, which is consistent with earlier suggestions 
of such pheno mena by | Dqnati et alJ J200ll) and the theoretical 
predictions of Gagn e et alJ J2005I) 

6 l Ori C holds a special place in our understanding of mag- 
netism in intermediate- and high-mass stars, because it is one 
of only two O-type stars in which a magnetic field has been 
detected and characterized at multiple epochs. It is also one of 
the youngest stars (with an age of about 10 6 years) in which 
a magnetic field has been detected, demonstrating once again 
(e.g. Bagnulo et al. 2004) that magnetic fields are apparent at 
the surfaces of some intermediate and high mass stars at very 
early main sequence evolutionary stages. The confirmation of a 
field in 9 l Ori C extends the range of known stars hosting fossil 
magnetic fields by a factor of about 2 in effective temperature, 
and by about 4 in stellar mass. 

Acknowledgements. GAW warmly acknowledges David Bohlender 
(Herzberg Institute of Astrophysics, Canada) for first bringing the in- 
triguing case of 8 [ Ori C to his attention. We also extend thanks to 
Otmar Stahl for helpful advice and fruitful discussions. GAW and 
IDL acknowledge Discovery Grant support from Natural Sciences and 
Engineering Research Council of Canada. 



References 

Babel, J. & Montmerle, T. 1997a, ApJ, 485, L29 

Babel, J. & Montmerle, T. 1997b, A&A, 323, 121 

Bagnulo, S., Hensberge, H., Landstreet, J. D., Szeifert, T., & 

Wade, G. A. 2004, A&A, 416, 1 149 
Baudrand, J. & Bohm, T. 1992, A&A, 259, 711 
Bevington, P. R. 1969, Data reduction and error analysis for the 

physical sciences (New York: McGraw-Hill, 1969) 
Charbonneau P. & MacGregor K., 2001, ApJ 559, 1094 
Cranmer, S. R., & Owocki, S. P. 1996, ApJ, 462, 469 
Donati, J.-E, Babel, J., Harries, T. J., et al. 2002, MNRAS, 333, 

55 

Donati, J.-E, Catala, C., Wade, G. A., et al. 1999, A&AS, 134, 
149 



Donati, J.-E, Howarth, I. D., Bouret, J.-C, et al. 2005, MNRAS, 
in press 

Donati, J.-E, Semel, M., Carter, B. D., Rees, D. E., & Collier 

Cameron, A. 1997, MNRAS, 291, 658 
Donati, J.-E & Wade, G. A. 1999, A&A, 341, 216 
Donati, J.-E, Wade, G. A., Babel, J., et al. 2001, MNRAS, 326, 

1265 

Fullerton, A. W 2003, in ASP Conf. Ser. 305, Magnetic Fields 
in O, B, and A Stars: Origin and Connection to Pulsation, 
Rotation, and Mass Loss, eds. L.A. Balona, H. E Henrichs, 
& R. Medupe, 333 

Gagne, M., Caillault, J.-P, Stauffer, J. R., & Linsky, J. L. 1997, 
ApJ, 478, L87 

Gagne, M., Oksala, M. E., Cohen, D. H, et al. 2005, ApJ, 628, 
986 

Johnson, N. M., 2004, M.Sc. Thesis, Royal Military College of 
Canada 

Landstreet, J. D. 1988, ApJ, 326, 967 

Mathys, G. 2001, in ASP Conf. Ser. 248, Magnetic Fields 

Across the Hertzsprung-Russell Diagram, ed. G. Mathys, 

S. K. Solanki, & D. T. Wickramasinghe, 267 
Moffat, A.F.J, and Michaud, G., 1981, ApJ 251, 133 
Reiners, A., Stahl, O., Wolf, B., Kaufer, A., & Rivinius, T. 

2000, A&A, 363, 585 
Semel M., 2003, A&A 401, 1 
Shore, S. N. & Brown, D. N. 1990, ApJ, 365, 665 
Shorlin, S. L. S., Wade, G. A., Donati, J.-E, et al. 2002, A&A, 

392, 637 

Simon-Diaz, S., Herrero, A., Esteban, C, & Najarro, F. 2005, 

A&A, in press (astro-ph/0510288! 
Smith, M. A. & Fullerton, A. W. 2005, PASP, 1 17, 13 
Stahl, O., Kaufer, A., Rivinius, T., et al. 1996, A&A, 312, 539 
Stahl, O., Kaufer, A., Wolf, B., et al. 1995, Journal of 

Astronomical Data, 1, 3 
Stahl, O., Wolf, B., Gang, T., et al. 1993, A&A, 274, L29 
ud-Doula, A. & Owocki, S. P. 2002, ApJ, 576, 413 
van Genderen, A. M., Alphenaar, P., van der Bij, M. D. P., et 

al. 1985, A&AS, 61, 213 
Wade, G. A., Bagnulo, S., Kochukhov, O., et al. 2001, A&A, 

374, 265 

Wade, G. A., Donati, J.-E, Landstreet, J. D., & Shorlin, S. L. S. 

2000, MNRAS, 313, 851 
Walborn, N. R. & Nichols, J. S. 1994, ApJ, 425, L29 



G.A. Wade et al.: The magnetic field and confined wind of the O star 1 Orionis C 

1 Ori C, 'x/v from longitudinal field variation, i=45° 
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Fig. 2. Map of reduced ^ 2 s of model fits to l Ori C longitudinal field variation (top) and to LSD profiles (bottom) assuming 
i = 45°. Reduced^ 2 intervals corresponding to 95% confidence [approximately 2cr) are computed using the^- 2 probability tables 
of Bevington(1969). 
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1 Ori C, LSD * 2 . i=25° 




Fig. 3. Map of reduced x 2 s of model fits to 9 l Ori C LSD profiles assuming i = 25° and i = 60°. Reduced^ 2 intervals 
sponding to 95% confidence are computed using the^ 2 probability tables of Bevington (1969). 
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Fig. 4. Averaged Stokes / and V and diagnostic null (AO LSD profiles of 6 1 Ori C, compared with profiles computed assuming 
the best-fit dipole model. 
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Fig. 5. Dynamic spectra illustrating the line profile variations of circumstellar features as a function of phase of the 15.422-day 
period. Quotient spectra computed by dividing each spectrum in the time series by the mean spectrum are illustrated. The original 
spectral time series is also overplotted in the lower panel, where the mean spectrum is shown in black. Note that different lines 
have different dynamic ranges ("stretches"). These are indicated by qmax and qmin, which are the maximum and minimum 
values of the quotient that are plotted, respectively. White (qmax) corresponds to wavelengths or times when the local flux is 
greater than its mean value. Strong nebular emission components have been excised from the central region of Ha and H/3 by 
interpolating linearly over -100 km s (for Ha) or ~50 km s _1 (for H/3). 
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Fig. 6. Same as Fig.|5J with a small dynamic range ("stretch') that is fixed for all features. This presentation emphasizes weak 
spectral features like the red-shifted emission near phase 0.5. 
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Fig. 7. Same as Fig.|5] only for a selection of features predominantly formed in the photosphere of 1 Ori C. 
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Fig. 8. Same as Fig. [5] only for a selected Hei lines. A strong nebular emission component has been excised from the central 
region of He i /15876. 



